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ignored, and will need to be included in future analyses. Here we provide the com- 
plete 0(p 4 ) expressions for these matrix elements obtained from Chiral Perturbation 
Theory, valid for partially quenched QCD with N degenerate sea quarks. Quenched 
QCD corresponds to the special case N — 0. We also discuss the role of the r{ meson 
in some detail, and we give numerical examples of the size of chiral logarithms. 
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1 Introduction 



The determination of non-leptonic kaon-decay amplitudes from Lattice QCD remains 
a challenging task. However, recently there have been several developments which 
may lead to substantial progress in this field, ranging from new ideas on how to cope 
with chiral symmetry on the lattice to a sizable increase of computer power available 
for the necessary numerical computations. 

Still, there are important theoretical difficulties afflicting the determination of 
the relevant weak matrix elements, which are a consequence of the fact that the 
final states contain more than one strongly interacting particle (the pions). This is 
formalized in what is sometimes called the Maiani-Testa theorem |l| , which says that 
it is not possible to extract the physical matrix elements with the correct kinematics 
from the asymptotic behavior in time of the euclidean correlation functions accessible 
to numerical computation. 

There are several old and new ideas on the market on how to deal with this situ- 
ation, which divide into three groups. First, one may determine the matrix elements 
with a kaon in the initial state and two (or more) pions in the final state for an un- 
physical choice of external momenta. From the large-time behavior of the euclidean 
correlation function 

C{t 2 ,h) = (0|7r 1 (t 2 )7r 2 (t 2 )C weak (t 1 )ir(0)|0) , (1.1) 

with the kaon at rest, one obtains the matrix element for p = p ni = —p- K2 = instead 
of the physical \p\ = (l/2)\Jm 2 K — Arri^, i.e. energy is not conserved fl], [|. The idea is 
then to use chiral perturbation theory (ChPT) in order to correct for this unphysical 
choice of momenta |J. The most recent computation of the AI = 3/2 K — > 2n 
matrix element using this method can be found in Ref. M. In this computation, 
all meson masses were taken degenerate and the quenched approximation was used. 
Adjustments for all these unphysical effects, which also include power-like finite- 
volume effects coming from pion rescattering diagrams, were made using one-loop 
ChPT [EJ. For recent ideas on choosing mx = 2m n for the lattice computation (for 
which p = does conserve energy), see Refs. || ||. 

For the AI = 1/2 case, the situation is more complicated for a number of reasons. 
Here we only mention (since this is less well known) that a quenched or partially 
quenched computation appears to be afflicted by "enhanced finite-volume effects," 
which do not occur for the AI = 3/2 case. This problem appears at one-loop in 
quenched ChPT f7|, |J, but not much is known about this effect beyond one loop. We 
are investigating this issue. For a nice review of many other issues, see Ref. [H]. 

A second idea was very recently proposed in Ref. ||, where it was shown how 
the matrix elements of interest can in principle be determined from finite-volume 
correlation functions without analytic continuation. The energy-conserving amplitude 
is obtained by tuning the spatial volume such that the first excited level of the two- 
pion final state has an energy equal to the kaon mass. With sufficient accuracy to 
determine the lowest excited levels, the finite-volume matrix element may then be 
computed on the lattice, and subsequently be converted into the physical infinite- 
volume amplitude. For this, it is obviously necessary to choose meson masses such 
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that 2m n < (as well as mx < 4:m n , so that the final-state pions are in the elastic 
regime). Again, if such a computation is done in a (partially) quenched setting one 
might expect that enhanced finite-volume effects could also occur with this method. 

A third idea is based on the observation that, if one needs ChPT anyway in order 
to convert an unphysical matrix element into a physical one, one might as well choose 
the unphysical matrix element as simple as possible. Chiral symmetry relates the 
K — ► 27T matrix elements of interest to the simpler K — > tt and K — > vacuum (K — > 0) 
matrix elements of the same weak operators which mediate non-leptonic kaon decay 
fT0|| . Advantages of this approach are that there are no strongly interacting particles 
in the final state, and that lattice computations of these simpler matrix elements may 
be less difficult. 

The first advantage is, in a sense, not really an advantage if one wishes to convert 
the results of a lattice computation into a calculation of the non-leptonic kaon decay 
rates, because final-state interactions will still have to be taken into account. However, 
this method does avoid all the unphysical effects, such as power-like or even enhanced 
finite- volume effects, associated with the multi-pion final state. Formulated in another 
way, with this method the simplest possible matrix elements (in this case K — > tt and 
K — > 0) are used to obtain the relevant weak low-energy constants (LECs) of the 
weak effective lagrangian. Using ChPT, these can then be converted into estimates 
of the kaon-decay rates. 

A key question is which order in ChPT will be needed in order to carry out 
such a program. At tree level {i.e. 0(p 2 )), only three LECs come into play, but at 
one loop (i.e. 0(p 4 )) many more LECs contribute to all relevant matrix elements 
Til . In fact, from the analysis reported in this paper as well as from previous work 



it is clear that tree-level ChPT is not enough 0, [7], |I2| . In addition (as we will 
demonstrate), not all 0(p 4 ) LECs needed for K — > 2tt decays can be obtained from 
K — > tt and K — > matrix elements. However, as we will advocate in this paper, it 
may be possible to determine at least the 0(p 2 ) LECs from a lattice computation, 
taking one-loop ChPT effects into account. A reliable, first-principle determination 
of the 0(p 2 ) octet and 27-plet LECs would clearly be interesting by itself. Moreover, 
phenomenological estimates of these LECs, based on a one-loop ChPT analysis of 
experimental data are available making a direct comparison possible. 

In this paper we present an analysis of K — > tt and K — > amplitudes in one- 
loop ChPT, with the above described philosophy in mind. For K — ► tt we choose 
our valence quark masses to be degenerate, thus conserving energy for this case. The 
analysis is performed in partially quenched ChPT ||14|| , with an arbitrary number of 



degenerate sea quarks. We present the results for these matrix elements in terms of 
the quark masses. 

In Sect. 2, we list and discuss all 0(p 2 ) and 0{p A ) operators needed for our 
calculation, including those containing the rf meson. In Sect. 3, we discuss the role 
of the rf in partially quenched QCD in some more detail than has been done so far 
in the literature. This section can be skipped if one is only interested in results. In 
Sect. 4, we give complete one-loop expressions for the octet and 27-plet K — > it and 
K —>■ matrix elements, including contributions from 0(p A ) operators, organized by 
subsection. In Subsect. 4.1 partially quenched results for N degenerate sea quarks 
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are presented, which are valid also in the case that the meson made out of sea quarks 
(the "sea meson") is not light compared to the rf (a realistic situation in actual 
lattice computations). In Subsect. 4.2 we specialize to the case that the sea meson 
is light compared to the vj '. Subsect. 4.3 contains the completely quenched results, 
obtained by setting N = and keeping the rf . For completeness, we include the fully 
unquenched results, with non-degenerate sea-quark masses for the K — > matrix 
elements, in Subsect. 4.4. In Sect. 5 we present a detailed discussion of the results, 
including the role of 0(p 4 ) operators and numerical examples for typical choices of 
the parameters. The last section contains our conclusions. 



2 Definition of operators 

Partially quenched QCD may be defined by separately introducing valence- and sea- 
quark fields, each with their own mass. The valence quarks are quenched by introduc- 
ing for each valence quark a "ghost" quark, which has the same mass and quantum 



numbers, but opposite statistics fllj, [15| . This, in effect, removes the valence-quark 
determinant from the QCD partition function. We will consider a theory with n 
quarks, of which N are sea quarks, and n — N > 3 valence quarks. This requires 
n — N ghost quarks, with masses equal to those of the valence quarks. We will 
consider valence quarks with arbitrary masses mi, . . . , m n _ N , and degenerate sea 
quarks, all with mass rris- The relevant chiral symmetry group is the graded group 



SU(n|n — N)L®S\J(n\n — N)r fl4| . Fully quenched QCD arises as a special case 
of this construction by taking N = fT6j, or equivalently, when the sea quarks are 
decoupled by taking m s — > oo. 

The euclidean low-energy effective lagrangian which mediates non-leptonic weak 
transitions with AS = 1 is given by 

Cas=i = C 2 + Ct + U + • • • , (2.1) 

where the dots denote higher order terms in the chiral expansion. The 0(p 2 ) la- 
grangian |l(J £2 contains three terms [] 

C 2 = - a * str (AL M L M ) + a\ str (AX+) + a 21 ' T^L^L^ 1 j + h.c. , (2.2) 

where str denotes the supertrace in flavor space. Note that the supertraces become 
normal traces, str — * tr , in the case N = n. The terms with couplings af, a\ 
transform as (8^, 1#), while the term with coupling a 21 transforms as (27^, The 
order p 2 lagrangian £?> w ^ be discussed toward the end of this section. The 0(p 4 ) 
lagrangian can be written as 

A = j^jy 2 (E A 8 Of + £ 61 + £ A 27 o? 7 ) , (2-3) 



the analysis of CP conserving weak amplitudes one can safely disregard O(e 2 p ) terms induced 
by electroweak interactions. 
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with (8 L , 1 R ) 0(p 4 ) operators Of, Of and (27 L , l R ) 0(p 4 ) operators Of JTTj, [T2J, |TT 
The operators Of denote total-derivative operators; they do not contribute to energy- 
momentum conserving matrix elements. However, they do contribute to the K — > 
matrix element, which does not conserve energy for non-degenerate quark masses. 
Note that there are no 27-plet total-derivative operators that contribute to the matrix 
elements considered in this paper, to 0(p 4 ). The fields entering the weak operators 
are defined as follows 

Lfi = zE^Et, (2.4) 
X± = 2S (EM t ± ME 1 ") , 



where B is the parameter B Q of Ref. |3i| and B = 4v/ f 2 in the notation of Ref. [|i~0 



The unitary field E is defined in terms of the hermitian field $ describing the Gold- 
stone meson multiplet as 

E = exp (2i$//) , (2.5) 

where / is the bare pion-decay constant, normalized such that f„ = 132 MeV. Note 
that L M and X± transform as (8l,1r) under (valence-flavor) SU(3)l®SU(3)#. The 
matrix A in the lagrangian (|2.2j ) picks out the AS* = 1, AD = — 1 part of the octet 
operators, all with AI = 1/2: 

A l , = 5% 2 , (2.6) 

where i, j = 1, 2, 3 . . . (or u,d,s...) denote valence flavors. The tensor Tjfi projects 
onto the AI = 1/2 part of the (27 l, 1#) operator, with nonzero components 

Til = ~ , (2.7) 



rpl3 

1 12 


rp31 

— 1 12 


rpl3 

— 1 2l 


rp23 

1 22 


rp32 

— 1 22 


= 1, 


rp33 

1 32 


rp33 

— 1 23 


3 

~ ~2 



or onto the AI = 3/2 part, with nonzero components 

rpV3 rp31 rpY3 T^31 ^ / o Q\ 

1 12 — 1 12 — 1 2l — 1 2l — r, J 
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rp23 _ rp32 _ _ _J_ 

1 22 — 1 22 ~ 2 



The term with coupling is known as the "weak mass term," and mediates the 
K — > transition at tree level. Its odd-parity part, which in principle can also 
contribute to the octet K — > tctt amplitude, is proportional to m s — m^. For m s ^ 
the weak mass term is a total derivative [0, [L9| , and therefore does not contribute to 
any energy- momentum-conserving physical matrix element, like K — > 7T7T. Instead, 
the K — > and, for Mk ^ M n , K —>■ n matrix elements do not conserve energy, and 
therefore the weak mass term does contribute to both of them. For m s = irtd the 
weak mass term is not a total derivative, so that it contributes also to the K — > n 
matrix element with Mk = M n . Hence, in order to determine the octet coupling 
af from a computation of the K — > tt matrix element, another quantity such as the 
K — > matrix element is always needed, in order to eliminate the dependence on a\. 
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At order p 4 there are eight (8^, 1 R ) operators and six (27 L , 1 R ) operators which 
can contribute to K — > 0, K — > 7r and K — > im matrix elements. The octet operators 



can be written as follows 



0\ = str (AX+X+) , (2.9) 

0\ = str (AX+) str (X + ) , 

Ol = str(AX_X^), 

Ol = str(A[X + ,X_]), 

O 8 10 = str(A{A + ,L /t L M }) , 

8 n = str (AL^X+L^) , 

8 n = str (AX+) str (L^LJ , 

0* 5 = str(A[A_,L^]), 



while the 27-plet operators are 



Of 


= tH(x + 


)\{x + ) l 3 , 


of 


= T%(X_ 


)\(X-) 1 J , 


of 


= T«(L, 


) i t ({L li ,X+}) l j 


of 


= TH(L, 




of 


= T%(X + 




of 


= nuL, 


)\(L,) 1 j str (X. 



(2.10) 



These operators are the same as those in Ref. |TjJ, apart from the replacement 



tr — > str . For the energy-momentum non-conserving matrix elements the only total- 
derivative term that is needed is 

6l = id fl str(A[L tl ,X + }). (2.11) 

In general, in the (partially) quenched formulation of the effective theory one 
needs to keep the rf QRj, |20], g| defined as the SV(n\n - N) L <g>SU (n\n - N) R 
invariant fl4j 

r/ = str($). (2.12) 



Note that this normalization differs from the one in Ref. |L4| , but is more convenient 



in keeping track of N dependence. The presence of the rf leads to new operators 



in the strong fig , [Tj| and weak effective lagrangians. Under CPS symmetry JT0[ all 
operators in £ 2 ,4 of Eq. ( |2.1| ) are even, and, since rf is CPS odd, new weak operators 
can be constructed by multiplying £ 2 4 by even powers of rf . It turns out that such 
operators do not contribute to the quantities of interest in this paper. However, other 
new weak operators arise from multiplying CPS-odd weak operators by odd powers 
of rf . There are two operators of interest at order p 2 , so that £ 2 is given by 

$ = 7i W//)str (AL„) + nl (t///)str (AX.) . (2.13) 

Since we are not interested in processes with external rf lines, we do not consider new 
operators of order p 4 containing the rf field. 



5 



For iV = 3 sea quarks the dynamics of the partially quenched theory is precisely 
that of unquenched QCD with degenerate quark masses JTJJ . Since all the low-energy 
constants (LECs) af, a 27 , 7® and /3j 8 ' 27 , pf (together with the strong counterterms) 
are independent of quark masses, it follows that their N = 3 partially-quenched values 



are equal to those of the real world ||23|| . However, for this equivalence to be valid, 
the rf should be treated in the same way in the partially quenched theory as in the 
real world, so that one has to consider the limit in which the rf decouples. 



3 The role of the rf in partially quenched ChPT 

Before we present our results, we would like to discuss the role of the rf in a partially 
quenched theory in more detail. First, define the bare (or tree-level) meson masses 

M 2 = B Q (m i + m j ) , i, j = 1, . . . 2n - N , (3.1) 

for a light pseudoscalar (pseudo-Goldstone) meson made out of quarks or ghost quarks 
i and j. For degenerate sea quarks, this simplifies to M| s = 2B ms. The two-point 
function for neutral mesons $jj (in that basis) is given by JE4] 



GM = J*rn*-xM> ( 3 - 2 ) 



where 

e» = 

and the rf mass is given by 



1 {m 2 + ap 2 ){p 2 + M 2 ss ) 



3 + Na {p 2 + Ml){p 2 + M?-)(p 2 + M 2 ) 

+ 1, for 1 < i < n 

-1, for n + 1 < % < 2n - 



]i = r + iv a /3 ■ (3 - 3 » 

The parameters and a (not to be confused with af' 27 ) come from the strong- 
lagrangian 0{p 2 ) operators quadratic in the rf field, (Nml/6)(r] r ) 2 + (Na/6)(d^rj') 2 . 
The term Xij(jp) has a double pole for = Mjj, unless M55 = Ma = Mjj. This 
implies that partially quenched theories suffer from the same "quenched infrared 
diseases" as the quenched theory unless all valence-quark masses are equal to the 
sea-quark mass [^4[ . From Eq. ( |3.2|) it is easily verified that the rf two-point function 
is just (N/(l + Na/3))(p 2 + M 2 )-\ 

The quantity Xy(p) can also be written as 



AiJ [P) a \p 2 + M 2 p 2 + M 2 + {p 2 + M|) {p 2 + Mfj)) [6A) 
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with 



M 2 - 



A; 



B; 



{N/?>){ml-aM%{m 



Ml s + ml{Ml s -M^{Ml s -M: 



331 



+ M? 



[(iV/3)K-«M| 

{N/3){m 2 -aM 2 ^ 



ss 
2 _ 




M|][(iV/3)K 



+ M 2 



ss 



+ a(Ml s -M? i )(M 2 ss 



+ M 2 



ss 



[(iV/3)K - <*M& + M 2 SS ~ MM(N/3)« 

{N/3){m 2 - aM 2 ss ) 2 /(l + aJV/3) 
[(N/3)(m 2 - aM 2 ^ + M 2 SS - M 2 } [(N/3)(m 2 - aM 2 ^ + M 2 SS - M 2 



3.P 



(3.5) 



The coefficients Ay, B^ and M 2 - are complicated functions of the various mass scales 
in the partially quenched effective theory. We may consider various limits in which 
these expressions simplify considerably. First, one easily obtains the fully quenched 
expression by setting N = 0, or equivalently taking Mss — ► 00, finding for all ij 



rn. 



■ 



.4 



1.1 



a 



B. 







(3.6) 



It is clear from these expressions that in the quenched case, the rj' should be kept in 
the effective theory. Another interesting limit is that in which the 7/ decouples p3] 
(which, as inspection of Eq. ( |3.3| ) tells us, is only possible for N > 0). In this limit, 
again for all ij, 



-Ms, 



A, 



(3.7) 



3 

N"~°° 1 ~~" J N ' 
|3.4j). The dependence of on the rj' parameters 



and we drop the 7/ pole in Eq 
has disappeared in this limit. 

As argued in Ref. [p0| , in actual partially quenched Lattice QCD computations 



the sea-meson mass Mss maybe comparable in size to mo so that the full dependence 
of A4 2 j and A^ on the parameters m , a and Mss should be kept. A third possibility 



1 j 



is then given by the limit in which the valence-meson mass is small compared to the 
7/ mass, i.e. M„ <C M„>. The expressions for Ai 2 jy Ai~ and B^ given in Eq. 



reduce to those of Ref. ]2(| if we expand in M% k /M 2 , but not in Mg S /M 2 ,: 



M 



1.1 



A; 



Br, 



M' 



A 



m 2 M 2 s 



(N/3)m 2 + M 2 S 



O 



'K, 



(3.8) 



{N/3)mt + aM, 



ss 



\{N/3)m 2 + M 2 s\ 



+ 



'Mi, 
Ml,, 



(N/3)(m 2 -aM 2 ss ) 2 /(l + aN/3) 
[(N/3)m 2 + M 2 s ] 2 



O 



'Mi, 
Ml, 



The (partially) quenched expansion we consider in this paper is systematic if we take 
A4 2 to be of order p 2 , in other words, if we take the parameter wig to be of the same 



order as the quark mass, just as in the case of quenched ChPT [[U| [2(J, |24 
It was also shown in Ref. 



20 



that, for the quantities considered there, simply 
ignoring one-loop contributions coming from the 77' pole in Eq. (|3.4|) and then taking 



the limit M v > 



00 in the rest is the same as matching to the limit in which the 
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Figure 1: A diagram contributing to the AS = 1 K —>■ n matrix element. A line is a 
pseudoscalar meson propagator, a dot a strong vertex and a box is a weak vertex. 



t]' decouples. In other words, if we ignore these contributions, the LECs appearing 
in those quantities are the same in the N = 3 partially quenched world and the 
real world. In addition, when we take only M^/M^ large, but not M^/Mss, these 
one-loop contributions are polynomial in the valence-meson masses, and can still be 
ignored if we are interested in the non-analytic dependence on the valence-meson 
masses. 

The same observations are also true here, even though there are diagrams with a 
more complicated topology than the simple tadpole diagrams needed in Ref. [p0|1 . For 



K — > 7r, there are contributions of the form depicted in Fig. 1. Taking M = M^k to 
be the degenerate valence-meson mass running in the loop and abbreviating B = Bkk, 
this diagram leads to one-loop integrals such as 



d 4 p 1 B B f d 4 p 



(2vr) 4 p 2 + M 2 p 2 + M 2 , M 2 , — M 2 J (2vr) 4 \p 2 + M 2 p 2 + M, 



(3.9) 



multiplied by two powers of M 2 from the two 0(p 2 ) vertices in the diagram (there 
are no contributions from vertices proportional to Mj 5 , so the only dependence on 
Mss comes through the coefficient B). The integral contributes an rf chiral loga- 
rithm, which we drop as above, and a Goldstone chiral logarithm proportional to 
BM 6 / (M 2 , — M 2 ) ■ log M 2 . If we do not assume that M is small compared to M v t this 
is of order M 4 , but if we expand in M 2 /M 2 ,, this constitutes an 0(p 6 ) contribution. 

In all the following calculations, we will assume that valence-meson masses are 
sufficiently small compared to M v i to justify the expansion in M 2 /M 2 ,. If one would 
not make this assumption, all contributions from integrals like Eq. (|3.9|) would have 
to be kept, since both the Goldstone-meson and the rj' pole give rise to additional 
non-polynomial M dependence. However, with this assumption, contributions com- 
ing from T]' tadpoles or diagrams such as Fig. 1 containing an 7/ on the loop are 
analytic in the valence-meson masses to order p A , and we will therefore drop them 
from consideration. 

Finally, we note that the coefficients of chiral logarithms will in general depend 
in a complicated non-polynomial way on the valence- and sea-meson masses, so that 
the 0(p A ) LECs cannot be defined in a mass- independent way. The coefficients of 
chiral logarithms have a polynomial dependence on meson masses only if we expand 
in both M 2 /M 2 , and M\ s jM 2 ,. In that case the LECs of the partially quenched theory 
without the rj' (and with N = 3) are the same as those of the real world. 
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4 K — > 7T and — > matrix elements at 0(p A ) 

In this section we calculate K° — > and K + — ► 7r + matrix elements at order p 4 in 
the effective theory. Four cases are considered: partially quenched with the valence- 
meson mass M <C M n > but M ss jM r] i arbitrary, partially quenched with also M v //M S s 
large, quenched, and unquenched. 



4.1 Partially quenched results 



We consider a partially quenched theory with three valence quarks, u, d and s and 
N degenerate sea quarks. The matrix elements to be calculated are defined as 



[K° 



7T 



0] 
+1 



(0|£ A 5=l|^°) , 
(n + \C AS=1 \K + ) . 



(4.1) 



The first process is calculated for m s ^ ma and m u = m^. In the last process we take 
the valence-quark masses all equal m s = = m u , so that it conserves energy. In 
this SU(3) limit, the K° — > 7] matrix element does not contain any extra information. 

We have calculated these matrix elements to 0(p 4 ) in partially quenched ChPT, 
using dimensional regularization in the MS scheme, and including contributions from 
the 0{p A ) operators ( |2.9| ) and ( |2.1U|) . In the case with degenerate valence quarks, let 
M = Mu be the physical mass of a meson made out of valence quarks, M$s that of a 
meson made out of sea quarks, and Mys that of a meson made out of a valence and 
a sea quark. At tree level in ChPT one has 



Mi 



vs 



;(M 2 + M 



ss) 



(4.2) 



Define, for any mass M, 



L(M) = log 



M 2 
A 2 



where A is the MS scale. In addition, define, for any two masses M\ and M2, 



L n (M u M 2 ) 



Mf n log ^ -M| n log 



Mi 

A 2 



Ml 



Ml 



(4.3) 



(4.4) 



As discussed in the previous section, we will assume that the valence-meson 
masses are small compared to M v / (cf. Eq. ( j3.3|) ), but not make the same assumption 
about Mss- At one loop, we then have, for the octet K — > n matrix element, using 
A4 2 and A from Eq. 



IK' 



7T 




(47T/)= 



NM VS 2 (L{M, 



vs 



(4.5) 



. 2{M 2 - t^AM 2 ) L(M) + ^(M 2 + 2AM 2 ) 



Aa\M 2 



(47T/) 5 



-(M 2 - AAM 2 )L(M) + -M 2 
3 3 
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4( 7l 8 + 2 72 8 )M 2 
/ 2 (4tt/) 2 



M 2 (2L(M) - 1) - ^NA (l 2 {M ss , M) - M 2 SS - M 
N M 2 - AM 2 



3 M 2 SS -M 2 



(2M 2 L(M) - L 2 (M SS , M) + M 



ss 



and for the 27-plet K — > 7r matrix element 
4a 27 M 2 / 1 



1" 27 



/ 2 



1- 



(47T/) 5 



2NM VS 2 (L(M VS ) 



1) + 2M 2 (3L(M) - 2) 



+ -(.M 2 - 2AM 2 )L{M) + -AM 2 
3 3 



(4.6) 



In the latter case, the matrix elements for AI = 1/2 and AI = 3/2 are the same, 
because of SU(3) symmetry. The r{ does not contribute directly to the 27-plet matrix 
elements; the dependence on m\ and a comes from the fact that we expressed all 
results in terms of bare meson masses, or, equivalently, in terms of quark masses (cf. 
Eq. (p.l|) ). The octet matrix element receives instead direct contributions from rj' 
exchange. Note that the contributions from the two 0(p 2 ) vj '-operators of Eq. ( p. 13 ) 
have the same form. This is explained by the fact that, after a partial integration, 
the first term in Eq. ( |2.13 ), using the equation of motion for E, is proportional to the 
second term. 

For the contributions from the 0(p 4 ) operators, we find 

' - ' \(2(3f + 2(3! + 2/4 + Pn)M 2 + (3 8 2 NM 2 SS (4.7) 



i(4) 



P 



(4tt/) 2 [^ 



+16a 8 ((A 4 - X G )NM 2 SS + (A 5 - A 8 )Af 2 ) - 8a 8 2 (X 4 NM 2 s + A 5 Af ; 
8M 2 - 



7T 



l(4) 
J 27 



P 



2(/? 2 + K ) M + ft NM ss 



(4.8) 



-16a 27 ((A 4 - X 6 )NM 2 SS + (A 5 - A 8 )M 2 



with again the AI =1/2 and AI = 3/2 results the same for the 27-plet. The Aj are 
the strong 0{p A ) LECs; they are related to the Gasser-Leutwyler Lj |18[ by 



A,- = 16rr L» 



(4.9) 



We see an example here of the fact that a partially quenched simulation, with 
M 7^ Mss, would in principle yield more information about the 0(p 4 ) LECs than 
an unquenched simulation in which M = Mss- 

For the K — > matrix elements we take non-degenerate valence quarks with 
m s 7^ m d and rrid 



m u , and define 



Ml, 



2M 2 K - Ml . 



(4.10) 



The pion will be made out of two light valence quarks, and the kaon out of a light and 
a strange valence quark. We also define M 2 S to be the (tree-level) mass of a meson 
made out of the z-th valence quark and quark, 



Mi 



B (mi + ms) 



u,d,s . 



(4.11) 
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Of course, M u s = M^s- For the octet matrix element we find, at one loop, 



[K - 0], 



N(M* s (L(M uS ) - 1) - M* s (L(M sS ) - 1)) 



(4.12) 



1 



+-M 1 M*(L(M n ) --)- M 2 3 (L(M 



33 / 



-A ( M*(L(M„) - \) - M| 3 (L(M 33 ) - | 



4mI(Ml - M 2 ; 



/ 



1 + 



(4tt/) s 



4 ( M us(HMuS 



M 2 sS {L{M sS ) 



2*7i 8 



--M 2 (L(M n ) + L(M 33 ) + 2L 1 (M 33 , M,) - 2) 
o 

+ (2M 2 L(ikQ + 2M 3 2 3 L(M 33 ) + 2L 2 (M 33 , M w ) - 3M 2 - 3M ; 

M w 4 (L(M 7r )-l)-M 3 4 3 (L(M 33 )-l) 

+^(M 2 - AM 2 ) (L 2 (M SS , M ff ) - M| s - M 2 ) 



iV 



--(.M 2 - AM 2 3 ) (L 2 (M S5 , M 33 ) - M 2 SS - M. 



4z 7 i(M 2 - Ml) 



2N 

Ml(L(M n ) - 1) + M 3 2 3 (L(M 33 ) - 1) - —AM 2 SS (L(M SS ) 



f{^f) 2 

+j(M 2 - AMDiL^Mss, M n ) - 1) + j(M 2 - AM 3 2 3 )(L 1 (M 5S , M 33 ) - 1) 

In this case, the leading order p 2 contribution comes only from the weak mass term 
proportional to a\. Unlike the case of K — > n, the contributions from the 0(p 2 ) 
rj'— operators proportional to 7^ and 7! are different in this case. The argument 
explaining the situation in the K — > tt case does not work here, because the total 
derivative in the partial integration cannot be dropped, as the process K — > does 
not conserve energy. 

For the 27-plet (both AI = 1/2 and AI = 3/2), we obtain 

12ia 27 T 

M*(L(M n ) - 1) + M 3 4 3 (L(M 33 ) - 1) - 2M 4 (L(M K ) - 1) 



27 



f(^f) 2 

2 



1 



+-M 2 MlL(M v ) + M 2 3 L(M 33 ) - L 2 (M 33 , M v ) + -(M 2 + M, 



33, 



~A (L 3 (M 33 , M w ) - ZMlM^L^M^ M n ) + 2M 2 M. 

Finally, the 0(p A ) operators of Eqs. ( [2.9| , [2.10| , |2.11|) give 

(4) 8t(M 2 K - M^ 



(4.13) 



f 



(2fi-2(3l + fi)M 2 K + (3»NM, 



ss 



-4a| (\ 4 NMl s + A 5 M 



12i 



(4nf) 2 [K^0}% = --f^l\M 2 K -Mlf 
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(4.14) 
(4.15) 



These results hold for an arbitrary number N of degenerate sea quarks. Note the 
appearance of "quenched chiral logs," contained in the one-loop logarithms propor- 
tional to M 2 . Since M 2 does not depend on the valence masses, such terms decrease 
with decreasing valence quark masses at the same rate as tree-level terms (modulo 
the logarithms), i.e. typically as mlogm instead of m 2 logm. 

Our results are presented here in a form somewhat different from that in Ref . . 
Here, we express the matrix elements in terms of the tree-level meson masses, or 
equivalently, the quark masses, while in Ref. |7j they were expressed in terms of 
renormalized masses (also, only the chiral logarithms were given). 

These results can be converted into expressions for the matrix elements as a 
function of the actual meson masses computed on the lattice by using the one-loop 
expression for the mass of a meson made out of two non-degenerate valence quarks 
in terms of the tree-level masses, Eq. Q3.1| ). This expression, in MS, and including 
0(p 4 ) contributions, is 



K) 



1— loop 



M 2 K [1 + 



(47T/) 5 



M 2 (L 1 {M 33 ,M 7r )-l) 



(4.16) 



+ -A(L 2 (M 33 ,M^) 



Mi 3 - m: 



+16 (M|(2A 8 - A 5 ) + NM 2 S (2X 6 - A 4 



For degenerate valence-quark masses, this simplifies to 



K 2 ) 



2\ 1— loop 



M 



1 



(4tt/) s 



--M 2 L(M) + -AM 2 (2L(M) - 1) 
3 3 

+16 (M 2 (2A 8 - A 5 ) + iVM| s (2A 6 - A 4 )) 



(4.17) 



4.2 Partially quenched results for large M„/ 

The results presented above simplify when M v i is taken large compared to both the 
sea- and valence-meson masses. Taking M v i large while keeping Mujj and M$s fixed 
in Eq. (|3.4|) gives 



1 



Mis 



M 2 - 

3 3 



N \p 2 + Ml {p 2 + Ml){p 2 + M 2 



(4.18) 



33 ■ 



dropping the 7]' pole. This expression for Xy(p) leads to the simplified expressions 
in Eq. (|3.7p for Ai 2 and A, where all the r/ parameters have disappeared, and the rf 
meson has been decoupled. Thus, for K — > n, we obtain 

Aa\M 2 ( 1 



[K- 



p 



Aa 8 M 2 



(4ir/) s 



NM VS 2 (L(M VS ) - 1) 



(4.19) 



N 



P 



((3M| S - 8M 2 )L(M) + M 2 S + 2M' 



(M 2 S - 4M 2 )L(M) + M; 



(Anf) 2 N 
12 



ss 



and 
[K+ 



4a 27 M 2 



71" 27 



P 



(4tt/) s 



2NM V s 2 (L(M\ 



vs) 



(4.20) 



+ 2M 2 (3L(M) - 2) + — ((M| 5 - 2M 2 )L(M) + M 2 



while for .K" — > 0, we find 

4iaf 



[K - 0]t 



iV(M 4 5 (L(M u5 ) - 1) - Mt s (L(M sS ) - 1) 



+-M 2 S (M 2 (L(M„) - -) - M 3 2 3 (L(M 33 ) - 
-| (m'MM^) -I)- M 3 4 3 (L(M 33 ) - |) 



(4.21) 



| 4iaf (M* - M 2 ) / + I 



./ 



(47T/)! 



N 



- (m 2 s (L(M uS ) - 1) + M 2 s (L(M s5 ) 



2iV 
1 



M| s (L(M n ) + L(M 33 ) + 2L 1 (M 33 , M v ) - 2) 



+ 



2A 

for the octet, and 

12ia 27 



(2M 2 L(M n ) + 2M 2 3 L(M 33 ) + 2L 2 (M 33 , M w ) - 3M 2 - 3M 2 3 ) 



27 



/(4vr/) 2 
2 



M^(L(M 7r )-l) + M 3 4 3 (L(M 33 ) 



2M 4 K (L(M K ) - 1) 



1 



+ j^M 2 ss (M 2 L(M W ) + M 2 3 L(M 33 ) - L 2 (M 33 , M^) + -(M 2 + M, 



~ (L 3 (M 33 , M w ) - 3M 2 M 33 Li(M 33 , Mjr) + 2M 2 M 



(4.22) 



for the 27-plet. All the dependence on rf parameters (m 2 ,, a, 7f 2 ) has disappeared 
from these expressions, as expected. The contributions from 0(p 4 ) operators, given 
in Eqs. ( |4.7|J4.8| ) and ( |4.14| , f4Tl~5| ), do not change. Note however, as mentioned before, 
that only in this limit the corresponding LECs are independent of the quark masses. 



4.3 Quenched results 

A special case of practical interest is the completely quenched result. In the quenched 
approximation, there are no sea quarks, and hence, quenched expressions can be 
obtained by setting iV = in Eqs. (|45 H4^j41^ , ^l^ , ^l^^T5|) , or equivalently, by 
taking M$s — * oo. In this case, it is not possible to decouple the rf |2T], [IE|| . The 
expressions given below maybe rewritten in terms of the parameter 5 (introduced in 
Ref. [H), by setting 

(4.23) 



24vr 2 / 2 5 . 
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The quenched results are 
Aa\M 2 l 



[K- 



1 



7T 8 



P 

Aa%M 2 



1 



(4tt/)2 

1 



/2 ^ ( 47r y)2 
4( 7l 8 + 2 7 f)M 2 



/W) 2 



2M 2 L(M) - M t 



2 (m 2 - gO;M 2 J L(M) + -(m 2 + 2aM 2 ) 
jj(m 2 - 4aM 2 )L(M) + jjm 2 ]) (4.24) 



4a 27 M 2 

T 2- 



1- 



(4tt/) s 



2M 2 (3L(M) - 2) 



2aM 2 )L{M) + -qM 5 



(4.25) 



[tf - 0] 8 



4zaf 
7(4^7? 



m 2 (M 2 (L(il4) 



i)-Af 3 2 3 (L(M33)-i; 



M 3 4 3 (L(M 33 ) 



4iaf (Mi - M 2 ) / 1 

/ I (4vr/) 2 

1 



(4.26) 



6 



m 2 Q {L{M n ) + L(M 33 ) + 2L 1 (M 33 , M ff ) - 2) 



+-a (2M 2 L(M 7r ) + 2M 2 3 L(M 33 ) + 2L 2 (M 33 , M n ) - 3M 2 - 3M 
o v 



2i 7l 8 



/(47T/) 2 

4i 7 f (M 2 - M 2 



/(47T/) 2 



M^ikQ - 1) - M 4 3 (L(M 33 ) - 1); 

M 2 (L(M.)-l) + M 3 2 3 (L(M 33 )-i; 



12ia 



27 



27 



/(4vr/) 2 



M*{L{M„) - 1) + M 3 4 3 (L(M 33 ) - 1) - 2M 4 K (L(M K ) 



+ jjm 2 (M 2 L(M 7r ) + M 3 2 3 L(M 33 ) - L 2 (M 33 , ikQ + ^(M 2 + M 3 2 3 ; 



-ia (l 3 (M 33 , Af ff ) + 3M 2 Mi 3 L 1 (M 33 , M n ) + 2M 2 M 



(4.27) 



Contributions from 0(p 4 ) operators are obtained by setting N = in Eqs. (4*77,4.8) 
and (|4.14| , [4.15|) . However, we emphasize that all the information obtained from 



quenched lattice computations is about the N = values of the LECs appearing 
in these equations. These values are, in principle, different from their N = 3 values. 



4.4 Unquenched results 

For completeness, we also report the results for the unquenched theory with three light 
flavors. For K — > tt these can be simply obtained by setting N = 3 and Mg S = M 2 in 
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our (large-M^/) partially quenched expressions of Subsec. 4.2. For K — > we have to 
choose the sea-quark masses equal to the non-degenerate valence-quark masses. The 
results for K — > therefore cannot be derived from our partially quenched results, 
where we took all sea quarks to be degenerate in mass from the start. The results are 

Aa\M 2 



71 



[K* 



7T 27 



P 

Aa 8 2 M 2 
4a 27 M 2 



1 + 



1 + 



(4tt/)2 

1 



(4tt/) 2 
1 

(477? 



-M L(M) + M 



2M 2 L(M) - -M 2 



-M 2 L(M) - —M 2 
3 V ; 3 



(4.28) 



for K —>■ TV, and 
[#->0] 8 = ^ 



(M 



A" 



Mi 



1 + 



~M 2 L(M.) - \m 2 k l{m k ) 



(4tt/)s 



(4.29) 



(4.30) 



U M 2 L(M V ) + -M 2 K + 



1* 



[jr-o] 



27 



/ (47T/) 2 

Aia 27 M 2 K - Ml 
~T (4vr/) 2 
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4z«? Mf M 2 \l_ L2{M ^ Mr) _ 2L2{M ^ + ^ + ^ M 2 



-2L 2 (M„, M K ) + 2L 2 (M ?? , Af w ) + 2M 2 K - 2M 2 



(4.31) 



for K — > 0. The contributions from 0(p 4 ) operators are obtained from the partially 
quenched expressions Eqs. ( fl~7jp||4TUlJ4.15|) , by replacing iVM| s -> 2M£ + M 2 in 
those equations. 



5 Relation to — ► 7T7T and numerical examples 

We now turn to a discussion on how our results can be used to extract physical 
information from lattice results. 

If tree-level ChPT were a good approximation, one could determine a\ 2 and a 27 
from a lattice computation of the K — > rr and K — > matrix elements, and then use 
ChPT to predict the AI = 1/2 and AI = 3/2 K — > rrrr decay rates [|Kj]. For instance, 
for the A/ =1/2 matrix element, one finds 

[K° - 7T + 7r-] 1/2 = 4 *(< -^) (a 8_ a 2 7) (51) 



<>K " " 



z mi- — m; 



[K + ^7T + ] 1/2 -6[K°^0 



/ M 2 

iM 2 _ 2im 
= f(M 2 K -M 2 ) = f(m s -m d ) ' ^ 

Here and are the physical kaon and pion masses, M is the degenerate meson 
mass (corresponding to a degenerate quark mass m) used in the lattice computation of 
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[K + — ► 7r + ]!/ 2 , and Mr and M„- are the non-degenerate meson masses (corresponding 
to quark masses m s and m^) used in the lattice computation of [K° — > 0]. At tree 
level, the procedure is very simple, because the conversion involves only meson masses 
and the meson decay constant /, which can relatively easily be determined. 

To repeat a similar analysis at one loop, one would not only need to eliminate 
the 0{p 2 ) constants a\ 2 an d a 27 , but also all the 0(p ) LECs that can appear in the 
matrix elements for the kaon decays of interest. In other words, the 0(p 4 ) weak LECs 



Pi 2 3 10 n 13 15 an d Pt 2 4 5 6 7 are needed |^9f| , as well as the strong LECs A^^s- Only a 
few linear combinations of those can be determined from K —>■ tt and K — > matrix 
elements. For the K — ► ttti matrix elements also only a few linear combinations 
are needed, but these are different linear combinations, involving also /3f 3 15 and (3f 7 



5,6 



which do not even appear in K — > tt and K — > at all. 

Also, we have seen that in unphysical matrix elements like K — ► new LECs 
appear in these linear combinations, such as for instance j3f in the linear combination 
j3f — 2/3f+/?f in Eq. (|4.14 ) . Likewise, one would be able to determine more LECs from 
K — > tt and K — > 77 in the mass non-degenerate case, but also more unphysical (total- 
derivative) 0(p 4 ) operators would contribute, since these matrix elements would also 
not conserve energy for onshell external states, just as [K° — > 0]. 

The only other way to determine more of the LECs from a lattice computation 
would be to consider more complicated correlation functions (such as K — > irir itself). 
In that case, one necessarily has more than one strongly interacting particle in the 
initial or final states, and this leads to rather severe complications of its own (for 
K — y niT, see Refs. [0, U). In general, on the lattice, one only has access to these 
matrix elements for unphysical choices of the kinematics 0, p], ||. Again, not all 
relevant 0(p A ) LECs can be determined. 

We conclude that at one loop uncertainties are introduced in the determination 
of K — > Till matrix elements from K — ► tt and K — > 0, which are not present at tree 
level. These uncertainties break down into two parts. One is the determination of 
a 27 from K — » rr and K — >• 0, and the other is the conversion of results for these 



a 



1,2 J 



0(p 2 ) LECs into K — > itit decay rates. Here, we will only consider the first part, i.e. 
we will concentrate on the determination of af 2 and a 27 from K —>■ tt and K —>■ 
matrix elements. 

If we are only interested in the determination of a\ 2 an d a 27 , we need to know 
only the polynomial form (in the meson masses) of the contributions from 0(p 4 ) 
operators, as given in Eqs. (|4.7[ , |4.8| , |4.14||4.15l) , assuming that one-loop ChPT can be 
reliably applied to lattice computations of K — >• tt and K ^ matrix elements. The 
results of the previous section can be written in the generic form 



[K -> tt] 8 



4M 2 



a 



I (l + Xf) -al(l + X%) + ( 7l 8 + 2 72 8 )X 

M 2 



(5.3) 



+Cv WJY 2+Cs WfY 2 



[K -> tt] 27 



4M 2 



a 



27 



(l + X 



_ 97 M 2 ~ N Mia 

m + cV—— + c 27 ss 



(4tt/)s 
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AT 
iV 


M 


X? 


X 2 


v27 
X 


3 


200 


0.72 


0.01 


0.34 


3 


350 


0.31 


-0.18 


0.74 


3 


500 


0.05 


-0.31 


1.12 


Z 


onn 
ZOO 


O.oy 


0.01 


n o a 
O.o4 


2 


350 


0.01 


-0.27 


0.74 


2 


500 


-0.47 


-0.47 


1.12 





200 


0.87(5/0.1) 


0.22(5/0.1) 


0.34 





350 


0.53(5/0.1) 


0.11(5/0.1) 


0.74 





500 


0.32(5/0.1) 


0.04(5/0.1) 


1.12 



Table 1: The one-loop corrections Xf 2 and X 27 for 2M|/3 = M 2 = M 2 , M ss = 
500 MeV, A = 1 GeV, M in MeV. The'iV = 2, 3 values have been calculated for large 
M v i, i.e. using the results of Subsect 4.2. 



[K - 0] t 



27 



Ai{M 2 K - Ml 
f 



Ai{Ml - Ml) 



f 



afY? + al (l + Yl 



M 2 

a 27 Y 27 + D 27 K 



M 2 



n s NMss 

+ s Wff 



(4tt/)= 



In these equations, Xf 2 , X 27 , Y® 2 , Y 27 , X 1 and Y\ 2 stand for the one-loop contribu- 
tions ("chiral logarithms") given explicitly in Eqs. (|4^ , |4.6|,|4.12lP^D . The constants 
Cy S etc. are linear combinations of 0(p 4 ) LECs, which can be expressed in terms of 



the /3's and A's by comparison with Eqs. ( |4.7|]4.8j , [4.14| , |4.15D . From these relations ct\ 



and «27 can be extracted by fitting these equations to lattice results for the matrix 
elements at various different values of the quark masses. With sufficient precision, 
also the 0(p 4 ) LECs could in principle be determined, but it is unlikely that this 
will work in practice with the currently available computational power. However, 
this does not imply that the 0(p 2 ) LECs cannot be determined with a reasonable 
accuracy. 

In order to get an idea about the size of one- loop effects, we will set the 0{p A ) 
coefficients Cyg , D vs and D 27 to zero, and evaluate the chiral logarithms at typical 
lattice values of the parameters and at A = 1 GeV, A = m p = 770 MeV, and 
A = = 550 MeV. We will consider three different "theories," partially quenched 
with N = 2 or 3 and M' n large, and quenched (N = 0) with arbitrary 5. We will also 
set f = fn = 132 MeV. We take Mss = 500 MeV, which corresponds to a sea-quark 
mass of about half the strange quark mass, vary the degenerate "lattice" meson mass 
M at which [K -> vr] is determined, and take 2M|-/3 = M 2 = M 2 for [K -»• 0], which 
corresponds to m s = 2m^ = 2m. 

It turns out that the one-loop correction X 27 for [K — > tt] 27 is very large. However, 
just as one defines the kaon B parameter Bk in the case of K° — K mixing, it makes 
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AT 
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0.00 


0.29 


3 


350 


0.23 
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500 
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n on 
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500 


-0.33 


-0.33 


0.83 





200 


0.71(5/0.1) 


0.17(5/0.1) 


0.29 





350 


0.37(5/0.1) 


0.06(5/0.1) 


0.60 





500 


0.16(5/0.1) 


-0.01(5/0.1) 


0.83 



Table 2: As in Table § but with A = 770 MeV. 
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Xf 
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X 27 
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0.42 
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500 
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500 


-0.14 


-0.14 


0.47 





200 


0.51(5/0.1) 


0.10(5/0.1) 


0.23 





350 


0.17(5/0.1) 


-0.01(5/0.1) 


0.42 





500 


-0.04(5/0.1) 


-0.08(5/0.1) 


0.47 



Table 3: As in Table 0, but with A = 550 MeV. 



sense to consider the ratio of this matrix element with its value evaluated by vacuum 
saturation or for large A^ c , which is proportional to (M 2 / 2 ) phys . In one-loop partially 



quenched ChPT, this can be expressed in terms of M and / as ||20 
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(m 2 / 2 ; 
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M 2 / 2 ll+A vs + 



>/) s 



(M 2 A 8 + NMlgXi 
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\(M 2 L(M) - AM 2 (2L(M) - 1)) 



+2NM VS 2 (L(M VS ) - 1) 



For the ratio 
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[K -> tt] 27 
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the relevant one-loop correction is 
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350 
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0.17(5/0.1) 


-0.16+0.07(5/0.1) 





500 


-0.10(5/0.1) 


0.10(5/0.1) 


-0.13+0.07(5/0.1) 



Table 4: The one-loop corrections F® 2 and Y 27 for 2M|/3 = M% = M 2 , M 55 = 
500 MeV, A = 1 GeV, M in MeV. The = 2, 3 values have been calculated for large 
M v >, i.e. using the results of Subsect 4.2. 



N 


M 


Yf 


Kf 


y27 


3 


200 


-0.63 


0.46 


0.01 


3 


350 


-0.44 


0.40 


-0.02 


3 


500 


-0.44 


0.41 


-0.02 


2 


200 


-0.42 


0.36 


0.05 


2 


350 


-0.07 


0.21 


0.02 


2 


500 


-0.02 


0.16 


-0.03 





200 


-0.36(5/0.1) 


0.23(5/0.1) 


-0.08+0.07(5/0.1) 





350 


-0.14(5/0.1) 


0.12(5/0.1) 


-0.09+0.07(5/0.1) 





500 


0.00(5/0.1) 


0.05(5/0.1) 


0.01+0.07(5/0.1) 



Table 5: As in Table [§ but with A = 770 MeV. 



In our examples, we will always consider the quantity X instead of X . From 
Eqs. ( |4.6|) and (|5.4| ) we see that X 27 is independent of N, M$s and the rj' parameters. 

In Tables |I] to |6| and Table [5] we computed the one-loop corrections for various 
choices of the valence-meson mass M and three values of the MS scale A = 1 GeV, 
770 MeV and 550 MeV. For = 0, we have set the parameters a and 7^ 2 equal to 
zero, for simplicity. In Figs. 2, 4 and 6, we show the dependence of the three largest 
octet one-loop corrections, Xf and Yf 2 on M and M55, for N = 2 and large M n i. 
We do not show similar plots for Af and Y 27 because these one-loop corrections are 
typically much smaller. We also do not show a plot for X 27 , because it only depends 
on M and not on M$s, nor on any of the f]' parameters. 

These examples illustrate various points: 

• One-loop corrections can be substantial, and will have to be taken into account 
in order to obtain a reliable estimate for af 2 and a 27 from the lattice. It 
may happen that the 0(p A ) LECs have values such that the combined, scale 
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AT 

iv 


M 




AX8 
Y 2 


Y 


3 


200 


-0.38 


0.32 


0.02 


3 


350 


-0.22 


0.26 


0.01 




500 


-0.27 


0.27 


0.05 


Z 


200 


-0.24 


0.25 


0.05 


2 


350 


-0.01 


0.13 


0.02 


2 


500 


-0.14 


0.13 


-0.03 





200 


-0.23(5/0.1) 


0.17(5/0.1) 


-0.05+0.07(5/0.1) 





350 


0.00(5/0.1) 


0.05(5/0.1) 


0.00+0.07(5/0.1) 





500 


0.14(5/0.1) 


-0.02(5/0.1) 


0.19+0.07(5/0.1) 



Table 6: As in Table § but with A = 550 MeV. 



A^ 


M 


M 


M{M) 


Moo 


A 


A(M) 


^4<x> 


3 


200 


433 


434 


500 


0.56 


0.61 


1.0 


3 


350 


433 


445 


500 


0.56 


0.73 


1.0 


3 


500 


433 


500 


500 


0.56 


1.00 


1.0 


2 


200 


499 


501 


612 


0.66 


0.74 


1.5 


2 


350 


499 


518 


612 


0.66 


0.95 


1.5 


2 


500 


499 


612 


612 


0.66 


1.50 


1.5 



Table 7: Comparison of values for M. and A, calculated from Eqs. ( |3.8|) (A4, A), (|3.4| ) 
(M(M), A(M)) and (.M^, A^), for 5 = 0.18, a = and M S5 = 500 MeV. M 
and are in MeV. 

independent contribution of 0(p 4 ) non-analytic terms and counterterms are 
smaller than the non-analytic terms alone at a given scale A < 1 GeV, improving 
the convergence of ChPT. One would hope this to be the case, especially for 
X 27 , which is very large at the larger values of M and A, and for Y± at larger 
A. This issue can be investigated on the lattice. 

• The size of the one-loop corrections grows with increasing A. In particular, they 
are relatively small for A = rs 550 MeV. However, without further knowledge 
of 0(p 4 ) LECs, it is unnatural to use m v as the scale, because the rj is itself a 
Goldstone boson, with mass very close to the meson masses we are considering 
here. In the absence of information on 0(p 4 ) LECs, we believe that using 
higher values for A gives a better a priori estimate of the size of 0(p A ) effects. 
If, however, the values of 0(p A ) LECs turn out to be such that A = 550 MeV 
gives the better estimate, one-loop corrections would be reasonably small, and 
one-loop ChPT should be applicable in the computation of af 2 and a 27 from 
the lattice at realistic values of the quark masses. 

• The fact that, in a number of cases, a one- loop correction becomes larger (and, 
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AT 

iV 


M 


X r 




V7 


3 


200 


-0.08 


A A A 

0.09 


A 1 /I 

0.14 


3 


350 


-0.12 


a 1 r 

0.15 


A 01 

0.21 


3 


500 


-0.17 


A 1 A 

0.19 


0.27 


2 


onn 
200 


-0.10 


0.11 


0.1 ( 


2 


350 


-0.16 


0.18 


0.26 


2 


500 


-0.21 


0.24 


0.34 





200 


-0.11 


0.14 


0.16 





350 


-0.23 


0.29 


0.35 





500 


-0.34 


0.40 


0.52 



Table 8: The one- loop corrections X 7 and Yfa for 2M|-/3 = M 2 = M 2 , M ss = 
500 MeV, A = 1000 MeV, 5 = 0.18 for N = 2,3, M in MeV. The results do not 
depend on 5 for N = 0. 



in fact, diverges) with decreasing M is a consequence of the fact that we do not 
vary Mss at the same time, i.e. of (partial) quenching. These "enhanced" chiral 
logarithms can be seen very clearly in the figures in the region M < Mss, toward 
smaller M. Chiral logarithms are typically smaller near the line M = Mss, f° r 
fixed M + Mss- Enhanced chiral logarithms appear in all quantities except X 27 . 

So far, we have used the large- M n * results, i.e. those of Subsect. 4.2, for all our 
examples. It is interesting to compare the JV^O values obtained at a finite rj' 
mass with those obtained for large M v '. In order to do this we chose 5 = 0.18 
for the finite- M^ case. (In the real world 5 ~ 0.18; for quenched QCD 8 ~ 0.1 



with a large error [ 25 1 ) . We set the other vj- related parameters a and 7{ 2 equal 
to zero, simply because not much is known about their values (but see below 
for some remarks on contributions proportional to these parameters). 

For N = 2, we show the difference between large and finite M v i in Figs. 3, 5 and 
7, for Xf and Yf 2 . We show the ratios [1 + X(finite M v )]/[1 + X(large M v >)], 
with X = Xf, Kf , and the difference (finite M v >) -Yf (large M^) (there is no 
tree- level contribution proportional to af). The first thing to be noticed is that 
there are still enhanced chiral logarithms in the small M, large Mss region. 
This is because the coefficients of the enhanced chiral logarithms depend on 

Mr,r. 

In the small M, Mss region, the ratios are close to one (or, the difference is 
close to zero), as one would expect if both M and Mss are small compared to 
Myf. For N = 2, M v > = 863 MeV, and M 2 SS /M 2 , w 0.34 for M ss = 500 MeV, 
while for N = 3, M v > = 996 MeV, so that M 2 S /M 2 , w 0.25 for M ss = 500 MeV. 
We point out that for larger meson masses, the plots are less meaningful in the 
region M > Mss- This is a consequence of the fact that we expanded the finite- 
M^i results in M 2 /M 2 f , but not in M^/M 2 ,. This makes sense if M < M ss , 
but not for M > Mss- Therefore, for valence-quark masses close to the sea- 
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quark mass, one either has to work consistently with the large- M v i results, or a 
more general analysis using Eq. (|3.4) instead of Eq. (p.8|) is needed. The latter 



is outside the scope of the present paper, but for Mg S /Mi ~ 0.25 — 0.34 it 
appears reasonable to use the large- M v i results of Subsect. 4.2. 

As a further illustration of this point, Table [7| shows the values of Ai 2 and 
A calculated from Eq. fl3T8f) (3rd, resp. 6th columns), the "exact" expression 
Eq. (|3.4|) with degenerate valence quark masses (4th and 7th columns), and 
in the limit of large M^, Eq. (|3.7|) (5th and 8th columns). We see that the 
variation of values for At 2 and A is at most about 20, resp. 50%. This is not 
very large from a practical point of view: the current best determination of 
Wq = At 2 (N = 0) in the quenched theory [^] has an error of the same order. 
We note that for M = 350 MeV the values of Ai and A are already much closer 
to their "exact" values Ai(M) and A(M) than .Moo and A^. 

All this means that, for the value of parameters chosen in our examples, one 
could use the simplest possible form of the chiral logarithms, given in Sub- 
sect. 4.2, to fit numerical results. Tables [1] to [5] then give the relevant examples 
of the size of non-analytic one-loop corrections. Only with numerical results so 
precise that one would be able to determine Ai 2 and A with better precision 
would it be important to take the dependence on rf parameters into account. 
Note again that these conclusions do depend on the values of the meson masses 
(and hence quark masses) we considered in our examples. 

We also considered the effect of the rf couplings jf 2 . They contribute to the 
octet matrix elements at one loop through the chiral logarithms X 1 and Y^ 2 . In 
Table [8] we show the size of these one-loop corrections for the choice A = 1 GeV. 
We see that these quantities are not very small, even though they vanish in the 
limit in which the rf decouples (for JV^O). They are typically smaller than Xf 2 
and Y^- If the couplings 7^ 2 themselves are small as well (compared to af 2 ), 
it may be possible to ignore these terms. But if these couplings are not small, 
this would mean that the rf does play a role, even if we can use the large- 
partially-quenched expressions of Subsect. 4.2 for the other chiral logarithms. 

The results are not very sensitive to a, except in the quenched case (N = 0). 
We expect that, within the precision of current lattice computations, the effect 
of (arbitrarily) setting a, as well as 7^ 2 , to zero, can be accommodated in fits 
by shifts in the 0{p A ) LECs C v etc. in Eq. (|5\3|), without significant impact 
on the values of the 0(p 2 ) LECs a\ 2 and a 27 . This can be checked by fitting 
lattice results while constraining these parameters to a few different values. 
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M [MeV] 



Figure 2: Xf (large M v >) with A = 1 GeV, as a function of M and M ss in MeV. 

6 Conclusion 

We presented a complete analysis of K — > 7r and — > weak matrix elements 
in one-loop ChPT for partially quenched QCD with N degenerate sea quarks. For 
K — > 7r we took the valence quarks degenerate in mass, while for K — > they are 
kept non-degenerate in order to get a non-trivial result. 

Three cases have been considered. The first is a partially quenched theory with 
a valence-meson mass much smaller than the r/ mass, but with arbitrary value of 
Mss/Mrf, so that an expansion in powers of M| fc /M2 is allowed. The second choice 
corresponds to the large limit where the rj' decouples, so that we can also expand 
in M^ s /M 2 ,. The third case is the quenched theory obtained for N = or Mss °o. 
For completeness, we also included results for the unquenched theory. 

These results should be useful for extracting the values of the LECs af 2 and 
a 27 from Lattice QCD. As we emphasized already in the Introduction, these are 
interesting quantities in their own right. Estimates extracted from experiment exist, 
with which lattice results can be compared. The matrix elements considered here are 
the simplest weak matrix elements which can be used for this goal. The expressions 
to be used in fits to lattice data are given in Eq. Q5.3Q , in which Xf 2 , X 27 , Yf 2 , Y 27 , 
X 1 and represent one-loop corrections (chiral logarithms), and Cy S , Dy S and 
D 27 are linear combinations of 0(p A ) LECs. They can be read off from the explicit 
one-loop results in Sects. (4.1-4). 

From our numerical examples in Sect. 5 it is clear that one- loop expressions will 
be needed for typical values of quark masses used in present lattice computations. 
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Figure 3: Rf = (1 + X* (finite M v ,))/(1 + Xf (large M v >)) with A = 1 GeV, as a 
function of M and M$s in MeV. 




Figure 4: (large M v >) with A = 1 GeV, as a function of M and M S s in MeV. 
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Figure 5: K x 8 (finite M v >) - (large M v> ) with A = 1 GeV, as a function of M and 
M ss in MeV. 




Figure 6: y 2 8 (l ar g e ^rf) with A = 1 GeV, as a function of M and Mss in MeV. 
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Figure 7: R\ = (l + y 2 8 (fmite M^))/(l+y 2 8 (large M n ,)) with A = 1 GeV, as a function 
of M and M ss in MeV. 

Contributions from 0(p 4 ) operators, represented by the LECs Cy 2 g , D vs and D 27 , 
will also need to be included. 

Theoretically, values for these 0(p 4 ) LECs can also be extracted from lattice 
computations. However, realistically, we expect that such estimates would have large 
uncertainties, both as a consequence of the typical statistics of present lattice com- 
putations, as well as because of uncertainties related to the role of the 7/ discussed in 
more detail in Sect. 5. We note that the 0(p 4 ) LEC /3f can be accessed directly by a 
computation of the ratio of the [K — > 0] 27 and the K° — K -mixing matrix elements 
(cf. Eq. (H)). 

In practice, for reasonably small values of the sea-quark mass (of order less 
than one-half times the strange-quark mass), it may be possible to use the partially 
quenched results for large M v i, given in Subsect. 4.2. In this case the rj' decouples, 
and therefore all dependence on rj parameters, m , a and 7 s 2 is removed, making 
the analysis simpler. In addition, it is only in this limit that estimates obtained for 
0(p 4 ) LECs can be directly compared to those of the real world, provided that one 
chooses N = 3 sea quarks. For a more detailed discussion, see Sects. 3 and 5. The 
more general results for the case that the sea-quark mass is larger, comparable to the 
t]' mass, but the valence-quark mass is still small enough, are given in Subsect. 4.1. 

A completely quenched lattice computation (for which the relevant results are in 
Subsect. 4.3) should be useful for assessing the feasibility of this approach. An N = 2 
computation, in combination with a quenched computation, could give insight into 
the dependence of the LECs on the number of light flavors. However, since we do 
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not know the functional form of the N dependence of the (finite part of the) LECs, 
an N = 3 computation will be needed to obtain estimates without an uncontrolled 
systematic error. 

The emphasis of this paper is on the extraction of reliable numbers for af and 
from lattice computations. While these 0(p 2 ) LECs are interesting, because of 



a 
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the availability of phenomenological estimates, the final aim of such Lattice QCD 
computations would be to convert these numbers into quantitative estimates of the 
AI =1/2 and AI = 3/2 K —>■ 2n decay amplitudes. This can be done using ChPT, 
and complete 0(p A ) formulae for doing so are given in Ref. ||19|| . Assuming that 0(p A ) 
ChPT is precise enough, the largest uncertainty arises because of the poor knowledge 
of all needed 0(p 4 ) LECs. Many of these, as we discussed in Sect. 5, cannot even 
in principle be determined from the K — > 7r and K — > matrix elements. (More 
0(p 4 ) LECs are accessible through the K° — K° and K — ► 2tt matrix elements with 
both pions at rest |J, the computation of which can also serve as a check on lattice 
results for af and a 27 .) One would have to resort either to the use of available 
phenomenological information [TBI, or to theoretical estimates based on arguments 



such as large N c or models (for recent discussions see Refs. jl^, |2l 

In addition, it is well known that final-state interactions are responsible for a 
large enhancement of the I = K — > tttt amplitude [|13|, |27|, p8[ . In that case it could 
be necessary to resum those effects instead of relying on an 0(p A ) ChPT calculation. 
A possible way of resumming final-state interactions has recently been proposed in 
Ref. IE7I. 
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